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ABSTRACT 



Context. Twisting motions of different nature are observed in several layers of the solar atmosphere. Chromospheric sunspot whorls 

and rotation of sunspots or even higher up in the lower corona sigmoids are examples of the large scale twisted topology of many 

solar features. Nevertheless, their occurrence at large scale in the quiet photosphere has not been investigated. 

Aims. The present study reveals the existence of vortex flows located at the supergranular junctions of the quiet Sun. 

Methods. We use a 1-hour and a 5-hour time series of the granulation in Blue continuum and G-band images from FG/SOT to derive 

the photospheric flows. A feature tracking technique called Balltracking is performed to track the granules and reveal the underlying 

flow fields. 

Results. In both time series we identify long-lasting vortex flow located at supergranular junctions. The first vortex flow lasts at least 
1 hour and is ~20" wide (~15.5 Mm). The second vortex flow lasts more than 2 hours and is ~ 27" wide (~21 Mm). 
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1. Introduction 

. Many measurements of the super granular flows have been car- 
X/^ ried out in the era of MDI/Soho (Shine et al. '20001; iPotts et alj 
-sj- 2004; Potts & Diver 2008; Meunier et al. 2008). However, the 
■rj" time averaging and spatial smoothing necessary to derive the 
flows could not allow their analysis at large scale with both high 
dynamic and high resolution. The new capabilities of the Solar 
Optical Telescope (SOT) onboard Hinode allow systematic long 
OO ^nd high resolution observations of the photospheric granula- 
tion. It is now possible to reveal more of the complexity of the 
L| underlying flows such as phot o spher ic vortex flows which were 
. ^ first observed by Bran dt et"aLl ( Il988l) . Our observations are de- 
scribed in section |2] followed in section [3] by a summary of the 
?H tracking algorithm used to derive the flows. Section |4] presents 
5^ , the results of this analysis. Finally we discuss some scientific 
implications of these observations. 



2. Observations 

We used two time series of de-rotated FG/SOT (Hinode) contin- 
uum images of the quiet Sun. They both have a field of view of 
~110"xllO"(~86x86 Mm^) : 

- 1-hour time series of Blue Continuum and Ca ii images on 10 
April 2007 between 17:00 and 18:00 UTC with a cadence of 
2 min at disk center. 

- 4-hours of contiguous G-Band and Blue continuum images 
on 07 Nov 2007 between 1:00 and 5:00 UTC near disk cen- 
ter at 3 min cadence. Ca ii images were also taken during 
the Blue continuum time series at the same cadence. We use 
also the Na I stokes V of the Narrow Band Filter in FG/SOT 
for preliminary estimates of the line-of-sight topology of the 
magnetic field. 



Both datasets were first binned onboard to 0.1 "/px (bin- 
ning 2x2) and corrected for dark currents, flat-fielding and in- 
strumental jitter using standard calibration procedures. For bet- 
ter computational efficiency, we binned the calibrated images by 
another factor 2 giving a final resolution of 0.2r7px and an im- 
age size of 512 X 512 px^. 



3. Analyses of photospheric flows 

To derive the photospheric flows from the continuum images, 
we used a t rackin g technique called Balltracking developed by 
iPotts et all (|2004|) . In this section we present a short summary 
of this technique. 

The balltracking technique considers the continuum images of 
the granulation in a three dimensional representation in which 
the intensity is seen as a geometrical height. Spherical tracers 
of a given mass and volume are released at a known position 
and settle in the local minima. In continuum images of the 
photospheric granulation these local minima are essentially con- 
stituted by the intergranular lanes. As the granules move, grow 
and shrink, the balls are pushed around and follow the resulting 
dynamic motion. The average width of the visible intergranular 
lanes is about 2 pixels (0.42") in our datasets which is the radius 
we set for the tracking balls. The outputs of the balltracking are 
velocity fields derived from the displacements of each tracking 
ball. They reflect the fast and stochastic small-scale granular 
motions. Therefore care is needed in the time averaging and 
spatial smoothing to reveal the underlying flows. Finally, a 
calibration procedure is needed to obtain physical units. As 
the balls take some time to settle in the local minima, the 
velocities are, on average, underestimated. The calibration 
consists in applying a known uniform velocity on subsets of the 
images to emulate an instrumental drift. The balltracking is then 
performed on the drifting subsets. Repeating the operation at 
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different value of the drift velocity reveals that the calculated 
velocity has a linear response and requires a correction factor 
of ~0.7. All smoothed balltracked velocities have been divided 
this factor. 

From simula ted velocity field and granular motions, 

iRieutord et al.l (l200lh observed that granular tracking gives 
relevant underlying flow fields only at spatial scales larger 
than 2.5 Mm and at temporal scales of at least 30 min. Thus 
observations of the velocity fields below these values may not 
be meaningful. At these thresholds. I Potts et al.l (|2004|) estimate 
that whatever the tracking algorithm used, for 1 min cadenced 
data and assuming perfect tracking of the granules, the uncer- 
tainty on the derived velocity flows would be ~150 m s ^ The 
uncertainty accounts for the finite scale of the granulation which 
makes the measurements spatially and temporally correlated. 
Thus the eff'ective number of independent data points is smaller 
than the number of samples taken by the instrument. Since our 
datasets have a cadence of 2 and 3 min, the uncertainty on the 
derived ve locities hav e been recalculated for each, following the 
formula in lPotts etal] (|2004 : 



cr,, = 



0"„ 



(1) 



where n,- and n, (respectively) are the efi'ective spatial and tempo- 
ral number of samples and cr„ is the standard deviation of a gaus- 
sian distribution approximatin g the velocity distrib utions of the 
granules. It is set to 800 m s"' dRoudier et al ] ll999h . Our veloc- 
ity fields were spatially smoothed over 4 Mm for both datasets. 
The results from the 1-hour and 5-hour observation series were 
averaged over 30 min and 45 min respectively. The choice of 
4 Mm, which is above the 2.5 Mm threshold, is a compromise 
between a better accuracy of the velocity measurements and a 
spatial resolution small enough to distinguish the motions in nar- 
row network lanes. Figure[T]represents the uncertainty on any of 
the 2 components of the horizontal velocity flow with respect to 
different instrumental cadences for a 45 min average. The uncer- 
tainty on the velocity flows averaged over 45 min from the 3 min 
cadence dataset is then ~78 m s '. At 2 min cadence and for a 
30 min average (not shown here), the uncertainty is very similar. 



- 


—2.5 Mm 
— 3 Mm 
-^d Mm 














- 
















- 




















- 
- 


















- 


















- 


















- 

































20 40 60 80 100 120 140 160 160 200 
Cadence (seconds) 

Fig. 1 : Uncertainty on the velocity measurements o",, for a 
45 minutes time averaging with respect to different instrumen- 
tal cadences. 



Nevertheless, considering that the mean lifetime of the granules 
is ~5 min, we believe that the granular motions cannot be prop- 
erly sampled at cadences of 2 and 3 min. Moreover the removal 
of the instrumental jitter needs at least 2 images during the time 
over which the features do not change significantly. It cannot be 



guaranteed in these datasets. For these reasons, quantitative re- 
sults may not be reliable. However, since we consider features 
of greater scale in time and space than the granular motions, the 
shape of the under lying velocity field can still be properly re- 
vealed jPotts et al.l2004,) . As a test, we used a 3-hour-time-series 
with a cadence of 20 s from which 3 datasets with a cadence of 
1,2 and 3 min were extracted. They were balltracked, smoothed 
and averaged over 4 Mm and 30 min. The upper pannel of figure 
12] shows that they are very well fitted by a Rayleigh distribution 
with a modal velocity decreasing with the cadence. Therefore 
we believe that lowering the cadence until 3 min acts as a low- 
frequency filter in the balltracking, which only picks up slower 
granules with lifetimes greater than 5 min. The bottom panel in 
figure |2]represents the distribution of the velocities of the April 
dataset. The modal and mean velocities are between 200 and 
300 m s"' . Th ese results are very c lose to the photospheric flows 
measured by Krii ger et al.l (12002 :). It is important to note that 
these values are smaller than the usual average speed measured 
in granular flows. A perfect tracking of the granular flow needs 
smaller smoothing and averaging. When we reduce these values, 
the modal and mean velocities indeed inc rease toward highe r 
values more consistent with the resul t s in [ Wang et al.l (Il995h . 
iBerger et al.1 ( [l998.) and iRoudier et al.1 (Il999h . To track flie un- 
derlying velocity flow, the smoothing and averaging must be in- 
creased to reduce the uncertainty introduced by the stochastic 
nature of the granular motions. 



— 60 3 Rayteigh fit 

• 120 3 

- - 1^ s Ftaytsi^h fit 
1SDs 

1B0 3 Rayteigh fit 




400 a» 
Vctacitv jm^s) 



VQiDchy ctata 
RaylBigh fit 




400 500 600 
Velocity (m/s) 



Fig. 2: Top: Distributions of balltracked velocities at 3 cadences 
for test data on 11* November 2007. Bottom: Distribution of 
the balltracked velocities for the 10* April 2007 dataset. Time 
average is 30 min and spatial smoothing has a FWHM of 4 Mm. 
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4. Results 

The derived velocity flows from the balltracking on both datasets 
are shown in figure [3] and |4l The upper (April) and bottom 
(November) row contain velocity fields at two different times. 
The black arrows are chained massless floating corks with an 
arrow head showing the direction of the flow (figure|5]l. 

These objects are updated each time the corks are at a differ- 
ent position in the stream. The light blue lanes in the background 
are the emphasized borders of the supergranula r cells derived 
from an automated pattern recognition algorithm jPotts & Diveil 
12008 *). The representation is similar to what floating corks would 
show once they gather in the converging sites of the flows. A 
large scale converging vortex flow is clearly visible in each 
dataset. The first one is ~20"(~15.5 Mm) and the second one 
is larger, ~27"(~21 Mm). Both of them are located at super- 
granular junctions as indicated by the blue network lanes. There 
relative position in each field of view is very similar but is only 
a coincidence. However, these observations were made at simi- 
lar latitudes and longitudes (near disk center). How it affects the 
occurrence of these vortex flows may be investigated with more 
observations and a greater statistic. Figure |5] sketches how the 
chained corks are dragged towards the vortex flow. 

Figure |6] is the Eulerian representation of the November ve- 
locity field centered on the vortex flow. The colored background 
is the normalized angular velocity. It is represented here as a vi- 
sual indicator and should only be considered qualitatively (see 
section 121). We do not see the beginning of the vortex flow in the 
first dataset and it seems to continue longer than the 1 h time se- 
ries. In the 2nd dataset, the onset seems to be between 2:30 and 
3:10 and continue past 05:00, the end of the time series. 
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5. Discussion 

We demonstrate the occurrence of steady photospheric vortex 
flow located at supergranular junctions of the quiet Sun chro- 
mospheric network. Their influence extends to a radius of at 
least 7 Mm from the center of th e vortex. The first observations 
of vortex flow was reported by iBrandt et al.l (1 1 9881) and had a 
radius of 2.5 Mm. At granular scales, excess of vortic ity has 
also been observed at th e intersections of granular lanes jZirkeij 
119931; IWang et aljri995[) . Magnetic bright points following spi- 
rals has recently been reported by Bone t et al] (12008). Long- 
lived vortices are of primary interest because these are the ones 
in which mixed polarity magnetic fields can become entwined. 
In the quiet Sun, magnetic fields are swept by the supergranu- 
lar flows to the boundaries and along the lanes to the core of 
the vortex, provok ing flux cancellation and CME-like eruptions 
dlnnes et alj|2008h . Fi gures |2] and [8] illustrate the possible con- 
nection between the magnetic flux and the vortex flow. In figure 
[8] the black and white colors represent opposite polarities of the 
line of sight magnetic field. The positive (white) polarity near 
the image center coincides with the Ca ii emission at the center 
of the vortex in figure |7] In the time series, it is possible to see 
the negative polarity (black spot pointed by the yellow arrow) 
moving clockwise around the positive polarity. 
Future observations must be done at higher cadence in or- 
der to reduce the uncertainty and to derive quantitative param- 
eters such as the angular velocity and the magnetic helicity 
dWelsch & Longcop3l2003h . Observations from Hinode and in 
the near future from SDO could be used to compute the build-up 
of energy and relate it to activity in the chromosphere and the 
corona. 
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Fig. 5: Sketch of the Lagrangian's representation of a vortex 
flow. The concentric circles represent the iso-contours of the tan- 
gential velocity increasing toward the center of the flow. The 
segmented arrows represent massless chained corks. 
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Fig. 8: Na I stokes V (FG/SOT) from the November dataset 
at 3 different times centered on the vortex flows. The FOV of 
~20x20 Mm^ is designated by the red square in fig.|7] The yel- 
low arrow points at a unipolar (black) magnetized fluid element. 
The dashed cross is a graphical fixed reference to help visualize 
the relative motions. 



Fig. 6: Eulerian representation of the flow in the red square of 
fig- El The colored background is the normalized vertical curl. 
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Fig. 7: Zoomed-in velocity field of fig. |4] with a filled contour of 
the maximum Ca ii emission. The red square indicates the patch 
represented in fig. [8] 



